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Abstract 

We examine how the binding of light (A < 8) nuclei depends on possible variations of hadronic 
masses, including meson, nucleon, and nucleon-resonance masses. Small variations in hadronic 
masses may have occurred over time; the present results can help evaluate the consequences for 
big bang nucleosynthesis. Larger variations may be relevant to current attempts to extrapolate 
properties of nucleon-nucleon interactions from lattice QCD calculations. Results are presented 
as derivatives of the energy with respect to the different masses so they can be combined with 
different predictions of the hadronic mass-dependence on the underlying current-quark mass m q . 
As an example, we employ a particular set of relations obtained from a study of hadron masses 
and sigma terms based on Dyson-Schwinger equations and a Poincare-covariant Faddeev equation 
for confined quarks and diquarks. We find that nuclear binding decreases moderately rapidly as 
the quark mass increases, with the deuteron becoming unbound when the pion mass is increased 
by ~60% (corresponding to an increase in X q = m q /h.QCD of 2.5). In the other direction, the 
dineutron becomes bound if the pion mass is decreased by ~15% (corresponding to a reduction of 
X q by ~30%). If we interpret the disagreement between big bang nucleosynthesis calculations and 
measurements to be the result of variation in X q , we obtain an estimate 5X q /X q = X-(0.013±0.002) 
where K ~ 1 (the expected accuracy in K is about a factor of 2). The result is dominated by 7 Li 
data. 

PACS numbers: 98.80. Cq,26.35.+c,21.45.+v 
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I. INTRODUCTION 

There are several reasons to search for a variation of the fundamental "constants" of 
nature in space and time. First, the Universe is evolving; several phase transitions happened 
in the early Universe accompanied by dramatic changes in vacuum energy, fundamental 
masses, fundamental interactions (electromagnetic, weak, and strong) and properties of 
elementary particles (e.g., confinement of quarks). At later stages the equation of state of 
the Universe continued to evolve, from radiation domination (pressure p = e/3 where e is 
the energy density), to cold- matter domination (p <C e), and "recently", about 5 billion 
years ago, to dark-energy domination (jp « — e). In view of these dramatic changes it seems 
natural to check if there is any evolution in the values of the fundamental constants during 
this process. 

A second reason is that spatial variation can explain a fine-tuning of the fundamental 
constants which allows humans (and any life) to appear. Indeed, it is well-known that if the 
values of some fundamental constants (e.g., related to the strong interaction) would differ 
by even 1% we could not appear. If we assume that the fundamental constants vary in 
space, this problem of fine-tuning may be resolved in a most natural way: we appeared in 
the area of the Universe where the values of the fundamental constants are consistent with 
our existence. 

A third reason comes from theories unifying gravity and other interactions. Some theo- 
ries suggest the possibility of spatial and temporal variation of physical "constants" in the 



Universe (see, e.g. jl, 2, 3, 0, S, S Q, Isj])- Moreover, there exists a mechanism for making all 
coupling constants and masses of elementary particles both space- and time-dependent, and 
influenced by local circumstances. The variation of coupling constants can be non-monotonic 
(e.g., damped oscillations). 

We can only measure the variation of dimensionless parameters which do not depend on 
the units we use. In the Standard Model the two most important dimensionless parame- 
ters are the fine structure constant a = e 2 /he and the ratio of the electroweak unification 
scale (determined by the Higgs vacuum expectation value) to the quantum chromodynam- 
ics (QCD) scale Aqcd (defined as the position of the Landau pole in the logarithm for 
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the running strong coupling constant, a s (r) ~ const/ In {Kqcd^/^c) ) . The variation of the 
Higgs vacuum expectation value (VEV) leads to the variation of the electron mass m e and 
quark mass m q which are proportional to the Higgs VEV. The present work considers ef- 
fects produced by the variation of X q = m q / Kqcd where m q = (m u + m^)/2 is the average 
current- quark mass. 

Up to now a majority of publications about temporal and spatial variation of the fun- 
damental constants have considered effects of variation of a. However, the hypothetical 
unification of all interactions implies that the variation of the strong interaction parameter 
X q = m q j 'Aqcd rnay be larger than the variation of the electromagnetic a. For example, 
the grand unification theories (GUTs) discussed in Ref. [l| predict 

^~35^ (1) 
X q a 

The coefficient here is model-dependent but large values are generic for grand unification 
models in which variations come from high energy scales; they appear because the running 
strong-coupling constant and Higgs constants (related to mass) run faster than a. Indeed, 
the electroweak (i=l,2) and strong (i=3) inverse coupling constants have the following de- 
pendence on the scale v and normalization point v§. 

ar\v) = arVo) + hln{v/v ) (2) 

In the Standard Model, 27r6j =12) 3 = 41/10, —19/6, —7; the electromagnetic or 1 = (5/3)a{ 1 + 
a 2 x and the strong a s = a 3 . In GUTs all coupling constants are equal at the unification 
scale, aj(z/ ) = c<gut- We may consider two possibilities: if we assume that only acur 
varies, then Eq. ([2]) gives us the same shifts for all inverse couplings: 

Sa^ 1 = da^ 1 = Sa^ 1 = Sa^j. . (3) 

We see that the variation of the strong interaction constant a^{v) at low energy v 
is much larger than the variation of the electromagnetic constant a, since 5a^/a^ = 
(a 3 /ai :2 )5ai : 2 /ai, 2 and a 3 > ai i2 . 

The second possibility is the variation of the GUT scale (z//z/o in Eq. ([2])). This gives 

Sot^/h = 8oT 2 x lb 2 = Sa^/h (4) 

Note that in this case the variations have different signs since b\ and 62,3 have different 
signs. However, we expect an even larger enhancement of the variation of a 3 {5a 3 /a 3 = 
(6 3 a 3 /6i 52 ai i 2)5ai, 2 /ai i2 ). 



The variation of mj 'Aqcd can be estimated from the definition of Aqcd- The running of 
a s near the electroweak scale is given by 



b s ln(v/A QCD ) 



(5) 



Let us take v = m z where m z is the Z-boson mass. The variation of Eq. (jSj) and relations 
above give 

Aim i \ \ i A/i i m ) C 6a(m z ) ( ^ 



5(m z /A QCD ) 



(m z /A QC D) b s a s (m z ) a s {m z ) a{m z ) a{m z ) 
The value of the constant C here depends on the model used. However, the enhancement 

1/a ~ 100 should make the factor C/a large. Note that m z (as well as m e ^ q ) is proportional 

to the Higgs VEV. 

If this estimate is correct, the variation in vn e ^ q j Aqcd may be easier to detect than the 
variation in a. The cosmological variation of m q / Aqcd can be extracted from the big bang 
nucleosynthesis (BBN), quasar absorption spectra and Oklo natural nuclear reactor data 
For example, the factor of three disagreement between the calculations and measurements of 



the BBN abundance of 7 Li may, in principle, be explained 
the level of ~ 10~ 3 — 10 -2 [lcj (see also recent work 
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the variation of m q / Aqcd at 
13 Q. The claim of the variation 



of the fundamental constants based on the Oklo data in Ref. 



141 ] is not confirmed by recent 



studies 15] which give a stringent limit on the possible variation of the resonance in 150 Sm 
during the last two billion years. The search for the variation of m e j Aqcd using the quasar 
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18| 



absorption spectra gave a non-zero result in Ref. [3] but zero results in Refs 
The present time variation of m e ^/ Aqcd can be extracted from comparison of different 
atomic 19( and nuclear 20, |2lJ clocks. The review of the recent results on the variation of 



a, m, 



J Aqcd and ra q j Aq C d can be found in Ref. 



221. 



As we mentioned above, one can measure the variation of the dimensionless parameter 
X q = m q /A QC D- Note that m q w 4 MeV < A QC d ~ 220 MeV. As a result, nuclear 
parameters (nucleon mass, reaction cross-sections, etc.) are determined mainly by Aqcd- 
Therefore, in all calculations it is convenient to assume that Aqcd is constant and calculate 
the dependence on the small parameter m q . In other words, we measure all masses in 
units of Aqcd and will simply restore Aqcd i n the final results. For example, the relation 
between the variation of the proton mass and quark mass 5m p /m p = 0.0645m q /m q should 
be understood as SX p /X p = 0.064:5X q /X q where X p = m p j Aq C d and X q = m q j Aq C d- This 
approach was formulated in Ref. 



Predictions for the relations between the variations of quark and hadron masses may 



come from a variety of approaches, such as chiral perturbation theory |23|, [24], [25|, [26( or 
Dyson-Schwinger equations (DSE) 27]. Nuclear binding energies and spectra would vary as 
a consequence of the hadron mass variation, affecting a number of physical processes, such 
as BBN and the Oklo phenomenon. Large enough changes might alter the stability of some 
nuclei, e^., unbind the deuteron, bind the dineutron, or even make A=5,S nuclei particle- 



stable 



29]. In this work we estimate the changes in binding for A=2-8 nuclei that result from 



small changes in hadronic masses. We also evaluate the effect of much larger changes on the 
two-nucleon systems, which may be relevant to present attempts to extrapolate results from 
lattice QCD calculations |3(|. We first consider how changes in meson, nucleon, and nucleon- 
resonance masses would alter some representative nucleon-nucleon (NN) Hamiltonians. We 
then solve for the energy of the two-body systems exactly and calculate variational Monte 
Carlo (VMC) estimates for the larger nuclei with these forces. We report our results as 
derivatives of the energies with respect to the different hadron mass changes, so that the 
results can be utilized with different predictions of the coordinated changes between quark 
and hadron masses. Finally, we utilize the DSE predictions for hadron mass-dependence on 
the quark mass as an example to explore the effects on nuclear spectra and BBN. 



II. NUCLEAR HAMILTONIAN 



We examine several Hamiltonians of the form: 



i i<j i<j<k 



(7) 



Here Ki is a nonrelativistic kinetic energy operator, is a two-nucleon potential, and 
is a possible three-nucleon potential. We consider three different Argonne models for Vif 
Argonne vu (AV14) and Argonne t>28 (AV28) from 1984 31], and Argonne v\& (AV18) from 



1995 
1995 



33]. 



. In conjunction with AV18, we will use the Urbana model IX (UIX) from 



The AV14 and AV18 models are conventional NN potentials, while the AV28 model 
has additional explicit A (1232) degrees of freedom. The AV14 and AV28 were constructed 



together and fit to the same np phase-shift solution WI81 of Arndt and Roper 



341 ] so they 



are phase-equivalent. The intention was to use them in parallel many-body calculations to 
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study the effect of including explicit A's in the nuclear Hamiltonian. In practice, AV28 



has proven difficult to use, so beyond the two-nucleon system, only some triton 



351 ] and 



nuclear matter [36] calculations have been reported. However, AV14 has been widely used 



in a variety of few-body [37 



calculations. The AV18 is 



381 ] and dense nucleon matter 
an updated version of AV14 containing charge-independence-breaking (CIB) terms and a 
complete electromagnetic interaction. AV18 was fit directly to 4,301 pp and np data in the 
1993 Nijmegen partial- wave analysis [4^], and AV18 and AV18+UIX have become standard 
Hamiltonians for ab initio calculations of light nuclei 41] and dense matter j^ . 



All three of the Argonne potentials contain electromagnetic (EM) interaction, long-range 
one-pion-exchange (OPE), intermediate-range attraction, and short-range repulsion, written 
as a sum of operator components: 

% = ^(%) + EKM + «J(r«) + «§(r<,-)]Q& • (8) 
p 

The number of operators Ofj is 14, 28, or 18, as indicated by the AVxx designation. The dif- 
ferent operators are discussed in detail below; here we summarize the differences between the 
models. The AV14 and AV28 potentials both use an average nucleon mass, while AV18 keeps 
separate proton and neutron masses, which introduces a small charge-symmetry-breaking 
(CSB) term into Ki. The Ki for AV28 also depends explicitly on the A mass. In AV14 and 
AV28, t> 7 is just the Coulomb interaction between protons (with a form factor); in AV18 the 
magnetic moment interaction, vacuum polarization, and other small EM terms are added. 

The AW part of OPE is the same in AV14 and AV28, and an average pion mass is 
used. In AV18 there is a weaker OPE coupling strength, a slightly different form factor, 
and the small charge-dependent (CD) terms owing to the difference between neutral and 
charged pion masses are kept. The intermediate-range attraction is due primarily to two- 
pion-exchange (TPE) processes; in AV14 and AV18 this feature is approximated by using 
the square of the OPE tensor function T(m 7r r) as a phenomenological radial shape for v p T and 
adjusting strength parameters of the 14 or 18 operators to fit data. In AV28, 14 operators 
with explicit NN-NA— AA couplings are added to those of AV14; twelve of these have OPE 
range and two have intermediate and short range. These produce much of the intermediate- 
range attraction explicitly through coupled-channel effects. The coefficients of the first 14 
AW operators are refit to the data but are smaller in magnitude than in AV14. 

The short-range repulsion may be attributed to the exchange of heavier p and u mesons 



with suitable form factors for finite-size effects, but in all the Argonne models it is treated 
phenomenologically. AV14 and AV28 use a Woods-Saxon radial shape with 14 strength 
parameters while AV18 has a slightly more general shape with 26 parameters. 



A. Potentials 

The OPE potential between nucleons can be written as: 

v n( r ij) — fwNN XijTi ■ Tj + XijTij . (9) 

Here = 3r zi r Z j — 7$ • Tj is the CD isotensor operator which contributes when the difference 
between neutral and charged pion masses is retained, and 



X ij = ±(x°+2X±), (10) 
x tJ = ± (x° - X±) , (11) 



X % = — ^ m<? Pi™*)* ■ °3 + T{mr l3 )S t] ] . (12) 

The Sij = 3<7j • TijUj ■ fij — Oi • <jj is the usual tensor operator and Y{mr) and T(mr) are the 
normal Yukawa and tensor functions 

Y{mr) = —t( r ) , (13) 

T(mr) = (l + ^ + J^y) Y(mr)t(r) , (14) 

where \x = mc/h and a short-range form factor £(r) has been incorporated that makes both 
Y(mr) and T(mr) vanish linearly at the origin. In AV18, the X 0,± are calculated with 
explicit m n o and m,± masses and the scaling mass m s = m w ±. In AV14 and AV28 an 
average m T = |(to^-o + 2m 7r ±) is used, so that X^ vanishes, and the scaling mass m s = m^. 
The coupling /2 7Viv= o.081 in AV14 and AV28, and 0.075 in the more modern AV18. 
The intermediate-range and short-range terms in the potentials are given by 

v p j(r l3 ) = PT\m^r) , (15) 



v 



P( rij ) = (S p + Q p r + R p r 2 )W(r) , (16) 



where W(r) is a Woods-Saxon function. For AV14 and AV28, Qp=Rp=0, while in AV18, 
there are boundary conditions on v p such that no more than two of the three S p , Q p , and 



R p are independent and free to be fitted for any given operator. The associated operators 
Ofj include fourteen charge-independent (CI) operators that are common to all the models: 



Of= 1M = [1, a t ■ aj, Sa, L • S, L 2 , L 2 (a, • a,), (L ■ S) 2 ] ® [1, r t ■ Tj ] . (17) 

Here L is the relative orbital angular momentum and S is the total spin of the pair. The 
AV18 model has four additional CD and CSB terms: 

0£= 15 ' 18 = [1, a t ■ a,, ,%] ® T VJ , [r zl + r ZJ ) . (18) 

These latter terms are small, but important for fitting the differences between current pp 
and np scattering data. 

There are fourteen additional operators in AV28 that involve explicit A degrees of free- 
dom. The first two are 

Of= 15 = (a t • SjOfo • Tj) + (Si • a^Ti ■ Tj) + H.c. , (19) 
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Of= 1G = S^in ■ T^ + SgiTi ■ Tj) + H.c. , (20) 



where Sj (Tj) is the transition spin (isospin) operator for particle % that changes a spin 
(isospin) | state to a | state. The generalized tensor operator is S\j = 3<7j • r^Sj ■ — <Ji ■ Sj 
and H.c. denotes the Hermitian conjugate. These operators are part of a generalized OPE 
contribution 

f TTl \ ^ X r 

vlfinj) = U NN U NA — -mc 2 Y(m nj )0^ + T{mr l3 )0™ (21) 

V Trig / o 

that produces NN-NA transitions. The 7riVA vertex coupling is taken as f^NA = ZfnNN- 
Other operators provide OPE contributions to NN-AA, NA-AN, NA-NA, NA-AA 
and A A- A A transitions, with an additional coupling f n AA = ^firNN being introduced. 
Intermediate- and short-range interaction is also added in the diagonal NA-NA and 



AA-AA channels. The full set of operators is detailed in Ref. [31 ]. 

The Urbana model of three-nucleon interaction (TNI) is written as a sum of two-pion- 
exchange P-wave and remaining shorter-range phenomenological terms, 

V ijk = V^ P + Vg k . (22) 
The structure of the TPE P-wave term is expressed simply as 

= E ■ Tj, Tj ■ T k } + C^X™ X%] [ Ti ■ Tj, Tj ■ T k \) , (23) 



eye 
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where is the same spin-space function of Eq. (|T2l evaluated with the average pion mass, 
J2 C yc is a sum over the three cyclic exchanges of nucleons i, j, and k, and = ^A^. 
The Ap^ and the strength of the V { f k term are determined by fitting the binding energy of 



3 H 33| and the saturation density of nuclear matter {42]] in the presence of a given Vij] the 



parameters of UIX were selected in conjunction with AV18. 
B. Dependence on hadron masses 

We consider changes in the Hamiltonian that would be induced by small changes in 
the hadron masses. However, we do not consider chan ges in the coupling constants that 



might occur due to mass-dependent loop corrections 23l. l24j] . Variations of the nucleon and 



nucleon-resonance masses alter the kinetic energy operator: 



h 2 

Ki = -7^V 2 + K - m N )c 2 , (24) 



where rrii is the average nucleon mass = \{rn p + m n ) for AV14, but may be iba part of 
the time for AV28. For AVI 8 the kinetic energy operator has no tua component, but does 
include a small CSB piece: 

Ki = ^ V 2 - ( mn ~ m p \ Tz . V 2 (25) 
2m N 1 2m N \m n + m p J 1 

where = 2m p m n /(m p + m n ). There is also a small dependence on the nucleon mass in 
the magnetic moment parts of the EM interaction in AV18. 

Changes in mAf and rriA also affect the energy expectation value for AV28 through effects 
of the NN-NA-AA coupled channels. This can be visualized in a simple closure approxima- 
tion (see Eq. (33) of Ref. j3l| and the accompanying discussion) for the TPE diagrams where 
one or two intermediate A's are excited, propagate, and then de-excited by the transition 
potentials: 

V e ff w 2Vnn-na^, — —, zVna-nn (26) 

hii + {rriA — rn N ) 

+ Vatat-aa p 9 , — -tVaa-nn ■ 

tii + A[m,A — mN) 

The two mean-energy denominators E\ and E2 would vary, in part, as the kinetic energies 
of the intermediate states: 

Ei w + , (27) 

2m n 2m a 



with k an average intermediate momentum. This physical effect can be approximated in the 
AV14 and AV18 potentials by multiplying the intermediate-range strength coefficients I p of 
Eq. <HM by factors (1 + 5n) and (1 + 5a)- The above equations represent only some of the 
terms contributing to the and wa dependence introduced by the coupled channels, so 
we will fix the 5n and 5 a terms for AV14 by requiring that calculated two-body energies have 
the same mass dependence as AV28 without these factors. We will take the same factors as 
approximately correct modifications for AV18 also. 

The dependence on the pion mass (which vanishes in the chiral limit m q = 0) can be 
obtained from the original pion-exchange interaction which is proportional to the second 
derivative of the exponential potential, VjVjfexp (— fir)/r] — see, e.g. js]. As a result, the 
variation of the pion mass affects the strength of the OPE through the m 3 dependence of 
in Eq. ffT2]) and through the range of the Yukawa functions Y(mr) and T(mr), Eqs. (I13H14I) . 
Note that the scaling mass m s is not allowed to vary (this is only a coefficient which does 
not vanish in the chiral limit m q = 0). The dependence of OPE terms on pion mass is thus 
the same as given in Eq. (10) of Ref. 23|]. Neglecting the small effect of our short-range 
form factor £(r), the function m 3 Y(mr) has a constant volume integral. If m is increased, 
m 3 Y(mr) will be larger inside r = 2/fi and smaller outside. However, m 3 T(mr) will be 
smaller at all values of r > 0. 

The generalized OPE of AV28 that couples to intermediate NA and AA states effec- 
tively incorporates considerable multi-pion exchange effects. We can approximate this real 
physical effect in AV14 and AV18 by assuming the intermediate-range potential has the 
mass-dependence v j(r) oc [m^T {m^r)] 2 . The closure approximation of Eq. (1261) that justi- 
fies this connection applies most clearly to the first six static (L-independent) operators Ofj 
of Eq. (|T7|) because of the closed algebra of the spin (isospin) and transition-spin (isospin) 
operators. However, the L-dependent operators can also be affected through the coupled- 
channels, so we will also consider changing these terms in the intermediate-range potentials. 
We will also consider changes to the residual Vj(r) terms in AV28. Thus we will study vari- 
ations attributable to the pion mass in four stages with our Hamiltonians: (1) change in the 
OPE part of the potential v % (and v^ 1 , etc. for AV28) only; (2) change in both v n and the 
static TPE parts v f -1 ' 6 (TPE-s); (3) these changes plus the L-dependent TPE parts v^~ 7,1A 
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(TPE-L); and (4) these changes in Vij plus the TPE part of the three- nucleon potential V^ P 
(TNI). 

Representing the dependence on heavier meson masses m p ~ m u « my is much more 
problematic because the Argonne models do not have explicit heavy-meson exchange. How- 
ever, the effect of a short-range interaction is determined mostly by the volume integral 
of this interaction, e.g., the volume integral gives us the strength constant B if we want 
to approximate a short-range interaction by a zero-range interaction BS(r). Therefore, we 
may approximate the dependence on masses of heavier mesons by a change in the range 
parameters r and a of the short-range Woods-Saxon potential 

W(r) = ^ — . (29) 

1 + exp[(r - r )/a\ 

used in Eq. (fT6]) . To keep the same relative variation of the volume integral for the Woods- 
Saxon potential and the meson exchange potential exp(— myr)/r we change the parameters 
at the rate: 

5ro = 5a = 2 5m v 
r a 3 m v 

In this case, as my increases, the range of the Woods-Saxon decreases while the value at 
the origin remains constant and the volume integral varies as m v 2 . 

The repulsive core of the three-nucleon potential V^ k may be attributed to a combination 
of multi-pion and heavy-meson exchanges, and also some relativistic effects [43j. Owing to 
its complicated nature and phenomenological treatment and the fact that it gives a small 
contribution to energy expectation values, we do not attempt to determine its variation with 
changes in the hadronic masses. 

III. ENERGY CALCULATIONS 

We calculate the energies of the 1 S (np) virtual bound state and the deuteron for AV28, 
AV14, and AV18 by direct solution of the two-body equations. The energies of 3 H, 3,4 ' 5 He, 
6 ' 7 Li, and 7 ' 8 Be are calculated for AV14 and AV18+UIX usin g va riational Monte Carlo 



(VMC) methods. The VMC method is described in detail in Ref. [4lj and references therein. 
Briefly, we construct suitably parametrized trial wave functions \Ev and evaluate the upper 
bound to the ground-state energy, 

E v = SlpflfXl > Eq (31) 
11 



using Monte Carlo techniques for the mult i- dimensional integration. The parameters in \Ev 
are varied to minimize Ey, and the lowest value is taken as the approximate solution. We 
use a variational wave function of the form 



\*v) 



i + E US? 1 



i<j<k 



sU(i + u K 



i<j 



(32) 



The Uij and U™ 1 are noncommuting two- and three-nucleon correlation operators, reflecting 
the spatial and operator dependence of Vij and V^, and S is a symmetrization operator. 
The form of the antisymmetric Jastrow wave function ^fj depends on the nuclear state under 
investigation. For the s-shell nuclei the simple form 



\*j) = Ufc(r ij )\$ A (JMTT 3 )) 

i<j 



(33) 



is used. Here the f c ( central (spin-isospin independent) two-body correlation func- 

tions and $a is an antisymmetrized spin-isospin state, e.g., 



|$ 4 (0000)) = A\ tpJ.pTnJ.7i) 



(34) 



with A the antisymmetrization operator. The for p-shell nuclei is more complicated; 

details are given in Ref. 41J. The two-body correlation operator is a sum of spin, 
isospin, and tensor terms: 

= E «p(»v)Q& , (35) 

p=2,6 

where the Ofj are the static operators of Eq. (ITT]) . The central / c (r) and noncentral u p (r) 
pair correlation functions are generated by a set of six coupled, Schrodinger-like, differential 
equations which include the and a number of embedded variational parameters. These 
parameters are optimized in the energy minimization and then kept fixed as the hadronic 
masses are varied. The correlations are adjusted, however, because the altered interaction 
is used to regenerate the correlations in each case. 

The energies we obtain are shown in Table [J and compared to experiment. The VMC 
method is reasonably accurate for s-shell nuclei, giving energies ~2-3 % less bound than exact 



methods like Faddeev or Green's function Monte Carlo for a given Hamiltonian [4lJ. It is not 
as good for p-shell binding energies, but many other features such as density distributions and 
electromagnetic moments are in good agreement. We anticipate that small changes in the 
binding energies induced by small changes in hadron masses will be accurately tracked with 
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the VMC calculations. The comparison between AV14 and AV18+UIX models indicates the 
importance of including a three-nucleon interaction to approach the experimental energies. 

We evaluate the mass-dependence of the energies of the two- and multi-nucleon systems by 
changing the hadron masses one at a time, increasing the value by 0.1% and calculating 
the resulting energy, and then decreasing by 0.1% and repeating the energy calculation. The 
VMC calculations follow the same random walk stored from the initial energy evaluation to 
reduce the contribution of statistical noise. Results given in the following three tables are 
the dimensionless derivatives of the energy with respect to changes in the hadron masses 

AS{m H ) = 6E/E . (36) 
omH/mH 

These results can then be combined with any given model for how the different hadron masses 
are correlated with the underlying quark masses to give a total binding energy prediction: 

5m H (m q ) 



E(m q ) = E(m qo ] 



l + X>£(m H ) 



where m qo is the physical current-quark mass. A specific example will be given below. 



(37) 



A. Two-nucleon energies 

The two-nucleon results for AS (win) are given in Table HT1 for all three Hamiltonians. A 
simple approximate relation between changes in the deuteron binding energy Q and virtual 
bound state energy e v is given by [lol |: 

Ae„(m H ) _ _^/Q_ 
AQ{m H ) yfc " 

Despite the wide range of values for AS (ran) in Table [TT], this relation is valid within a 
factor of ~2 for fifteen out of seventeen pairs of results. The two exceptions are for the OPE 
components of Ae v (m n ) for AV14 and AV18, where the sign is incorrect. This discrepancy 
will be discussed below. 

Changing the nucleon mass in AV14 and AV18 primarily changes just the kinetic energy 
component and is very similar for the two models: a larger translates to smaller (K) 
and more binding and thus a positive derivative AQ(rriN) for the deuteron. For the virtual 
bound state, the effect is the same but greater attraction corresponds to a reduction in the 
magnitude of e v and a negative Ae v (rri]y). 
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Changing in AV28 produces a larger change in the energies due to the additional 
coupled-channel effects of Eqs. (I26H28P discussed above. The change is more dramatic for 
e v than for Q, probably because the 1 So NN channel can couple to both iVA and A A 
intermediate states, while the deuteron can have only NN and AA components due to isospin 
conservation. To approximately incorporate this physical effect into the phase-equivalent 
AV14 model, we can multiply its intermediate-range strength parameters I p of Eq. (fl5|) by 
a factor (1 + 5n)- We choose a value 5n = 0.495mjv/"Tjv that matches the mass dependence 
in the deuteron as shown by the line in Table [Til labeled + 8n" . This same factor 
approximately fixes the virtual bound state also. Using the same 8n factor in AV18 produces 
a change about one third larger. This difference is probably due to the fact that AV18 is 
fit to a more recent data set, with a weaker OPE coupling accompanied by more- attractive 
intermediate-range terms and correspondingly more-repulsive short-range terms. 

Increasing the A mass in AV28 reduces the binding energy both through the one-body 
term of Eq. (124]) and through the coupled-channel effects of Eq. (l2~6l) . Thus A£(7Ba) has a 
sign opposite to AS(rriN). To approximate this affect for AV14, we can again multiply its 
intermediate-range strength parameters I p by a factor (1+8 a)- We take 8a = — 0.57<5mA/mA 
which gives a fair reproduction of the behavior of AV28 for both the deuteron and the virtual 
bound state as shown by the line in Table UTI labeled "5a"- Using the same factor in AV18 
again produces a larger rate of change. 

We note that if the nucleon and A masses change at about the same rate, from what- 
ever the underlying quark mass dependence, then the effect on the AV28 energies could be 
obtained by the sum of the and wa terms in Table UTI which is not very different from 
the tun term alone for AV14 and AV18 or from the sum of + <5jv and 5a terms. In fact, 
these corresponding sums of terms agree within 10-20% for all three models. 

Increasing the pion mass reduces the regularized OPE tensor function rn^T(m^r) at all 
values of r. The binding of the deuteron is primarily due to the tensor coupling, so Q is 
reduced and AQ is negative for all the models as shown by the line in Table [TTJ labeled 
u m n (OPE)". The value for AV28 is larger because of the multi-pion-exchange effects in- 
cluded through the generalized OPE potential. The AV14 and AV18 values are smaller and 
close together, the difference between them being proportional to the different f% NN coupling 
constants used in the two models. 

Increasing m n in the 1 Sq channel, which depends only on m^Y (m^r) , has the effect of 
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slightly increasing the attraction in the virtual bound state, making A(e v ) negative for 
the AV14 and AV18 models. However, in the AV28 model, the generalized OPE provides 
significant intermediate-range attraction oc [m^T(m n r)} 2 , which is enough to reverse the 
sign of A(e t) ) to be positive. As the AV14 and AV18 cases here are the only two out of 
seventeen in Table HT1 that violate the relation between AQ and A(e„) of Eq. (1381) . it appears 
that changing only the OPE part of a conventional NN potential like AV14 or AV18 is an 
incomplete representation of the physics in the singlet channel. 

Changing m„ in the static TPE part of the interaction has a more dramatic effect than 
changing just the OPE part, as shown in Table HT1 by the line labeled u m n (+TPE-s)." The 
AQ increases by by a factor of 2.5-4 in the deuteron for the AV14 and AV18 models compared 
to OPE only, while for AV28 it increases by a smaller factor of 1.5. More dramatically, the 
AV14 and AV18 values for Ae v change sign and come into agreement with Eq. (138]) . The 
further addition of the non-static u m n (+TPE-L)" terms makes no difference to the virtual 
bound state, as these operators do not contribute in the 1 So channel, and rather small 
changes in the deuteron. 

Finally, increasing the heavy-meson mass my reduces the range of the Woods-Saxon 
repulsion and increases the binding of the deuteron, so AQ is positive. The reduced repulsion 
for the virtual bound state correspondingly makes Ae v negative. In both the full-pion and 
heavy-meson exchanges, the AV18 has larger AS values than AV14 or AV28. Again, this is 
probably due to the weaker OPE and compensatingly larger intermediate-range attraction 
and short-range repulsion. 



B. Multi-nucleon energies 

The multi-nucleon results for A£(m#) are given in Table [TTT1 for AV14 and Table [IV] for 
AV18+UIX. For every nucleus, and for every itlh component, the signs are the same as for 
the deuteron. The relative sizes of the terms are also about the same as the deuteron, with 
the exception of the m n (OPE) term. Because the light nuclei have approximately equal 



numbers of deuteron-like and 1 So-like pairs 44|, the A£(m#) are expected to be averages 



of the two sets of trends in Table HT1 The anomalous behavior of Ae v in the m n (OPE) case 
discussed above causes these numbers to be much smaller in the multi-nucleon systems, and 
thus out of proportion compared to all the other terms. As for the deuteron, the multi- 
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nucleoli terms have a more rapid dependence with AV18+UIX than with AV14. However, 
the explicit m n contribution from the three-nucleon force is very small. 

IV. DEPENDENCE OF NUCLEAR BINDING ENERGIES AND BIG BANG NU- 
CLEOSYNTHESIS ON QUARK MASS 

As an example of how to incorporate our nuclear binding energy results with a specific 
prediction for hadronic mass variation, we utilize the results of a Dyson- Schwinger equation 



study of sigma terms in light-quark hadrons [27[. Equations (85-86) of that work gives 
the rate of hadron mass variation as a function of the average light current-quark mass 
m q = (m u + m d )/2 as: 

rriH m H m q 

with UH/mtH values of 0.498 for the pion, 0.030 for the p-meson, 0.043 for the a;-meson, 
0.064 for the nucleon, and 0.041 for the A. The values for the p and cj-mesons were reduced 



to 0.021 and 0.034, respectively, in a subsequent study [28(. We will use an average of the 
p and <jj terms of 0.03 for the variation of our short-range mass parameter my. 

In Fig. [I] we show direct calculations for the two-nucleon states, i.e., the deuteron and 
1 So (rap) states, for a range of m q values for all three Hamiltonians. The virtual bound state 
energies are plotted as a positive quantity; when the energy is negative it indicates a 
true bound state. The dashed lines show the results from changing only the pion mass in the 
OPE part of the interaction, corresponding to the line m n (OPE) in Table HT1 To evaluate 
over this wide range of m q , we use the Gell-Mann-Oakes-Renner (GMOR) relation oc m q . 

The solid lines show the most complete calculation for each Hamiltonian. For the AV28 
model this is the sum of the terms m^, ffiA, m n (+TPE-L), and my, using GMOR for the 
pion and the DSE values above for the variation of all other mu- For the AV14 and AV18 
models, it is the sum of the terms m^ + 5n, £a, m^ (+TPE-L), and my. As discussed 
in Sec. III. A, the 5^ and 5 a effects are incorporated by multiplying the intermediate-range 
strength parameters I p of Eq. (|T5|) by a factor 

(1 + S N )(1 + 5 A ) = (1 + .49^)(1 - .57^) (40) 

m N m A 

which is ~ (1 + .008 5m q /m q ) for the DSE values above. 
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At the OPE level, the AV14 and AV18 models show almost exactly the same behavior for 
the 1 So state, with a gradually increasing attraction as m q increases, while for the deuteron 
they both show a more rapid decrease in binding, consistent with the results obtained in 



Refs. 
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24j. If these trends continue, the deuteron will eventually move above the sin- 



glet state somewhere in the range 3-4 m q . The AV28 model has a somewhat more rapid 
dependence for the deuteron but its singlet state parallels the deuteron, becoming less at- 
tractive for larger m q . In the other direction, the singlet state becomes a true bound state at 
~ 0.7m q , but it always remains above the deuteron. This different behavior is a consequence 
of the multi-pion exchange that is built into the AV28 model through the generalized OPE 
coupling to intermediate A's. 

In all these models the deuteron is bound largely through the tensor coupling between the 
3 Si and 3 Di AW states. In the AV28 model, the singlet state gets a considerable part of its 
attraction through the tensor coupling between the 1 So AW and intermediate 5 D iVA states, 
and thus has a sensitivity to changes in the pion mass similar to the deuteron. This behavior 
of the singlet state, i.e., that it parallels the deuteron, is different from that predicted by 
chiral perturbation theory evaluated at the next-to-leading-order (NLO) in Ref. 30(, which 
does not include the effect of the A. We expect that a higher-order chiral perturbation 
calculation that includes A degrees of freedom will come into qualitative agreement with 
our result. 

For our most complete calculations, shown by the solid lines, the mass dependence of the 
energies is significantly steeper for both deuteron and singlet states, but they all are parallel 
and it appears the deuteron will remain the ground state for a very large range of m q . 
The AV28 curves shift relatively little from the OPE-only values, while the AV14, which 
is essentially phase-equivalent with AV28, gives very similar results. The biggest change 
and the most-rapid dependence is given by the AVI 8 model, with the most important 
contribution coming from its static two-pion-exchange terms, i.e., the line m % (+TPE — s) 
in Table [III The more rapid dependence is a consequence of the deeper intermediate-range 
attraction and stronger short-range repulsion, which in turn may be a consequence of the 
improved quality of AV18, i.e., that it is a better fit to more recent AW data. 

The two-body energies can also be evaluated using Eq. (I3TI) to combine the DSE values 
for bmujmii and the dimensionless derivatives A£(m#) of Table HT1 This might not be 
expected to work for large changes in m q , or where a state is barely bound or unbound. 
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However, deuteron energies are reproduced for changes of ±0.1 in Sm q /m q to 1% or better 
by Eq. (J37D. 

The dependence of the multi-nucleon energies for the full calculation with DSE values is 
shown in Fig. [2] for AV14 and in Fig. [3] for AV18+UIX. These results have been calculated 
using Eq. (157)) and the A£(m#) of Tables HTT1 and HV) respectively. The lines have been 
extended to 5m q = ±0.2m q to show the trends, although the results are not expected to be 
completely linear over such a broad range. The values have been checked in a few cases by 
doing direct calculations and adjusting the variational parameters to reminimize the energy. 
For example, the 4 He energy can be lowered by 0.3-0.5 MeV at either end of its line, but 
this change would hardly be visible at the scale shown. 

The multi-nucleon energies parallel the deuteron, with generally increasing slope as the 



binding energy increases. (Similar results for the triton were found in Refs. 45|, |46j.) The 
curves are steeper for AV18+UIX than for AV14. In either case it appears that the relative 
stability of all the nuclei will be preserved across a broad range of m q values, with the 
exception of 8 Be. Here it appears that a decrease in m q of ~ 0.5% will lead to 8 Be stable 
against breakup into two a's for both Hamiltonians. A very weakly bound (~ 0.1 MeV) 8 Be 
might not have much of an impact on primordial nucleosynthesis because it would be easily 
photo-disintegrated until quite late in the BBN epoch. A moderate binding (~ 1 MeV) 
could have a significantly more dramatic effect on the chemical evolution of the Universe 
by giving rise to the production of noticeable amounts of stable elements with A=9,10. If 
the bound state persisted to the era of star formation, it would presumably also have a 
significant effect on stellar evolution. 

A summary of sensitivities of nuclear binding energies to the quark mass m q , as given 
by the DSE hadronic mass variation, is presented in Table M for the different Hamiltonians. 
The total sensitivity K, 

of deuterium binding energy to the light quark mass is = —1.39 for the AV18 interaction, 
while the pion contribution ranges from KJ = —0.70 for the OPE contribution only to 
K% = —3.36 when the full TPE is counted. This may be compared with the pion contribution 



3 > K$ > -18 from Ref. 
The result of Ref. [l0| 



9] and = -2.4 from Ref. [24 1. 
suggested that a reduced deuteron binding energy of AQ = 
—0.019 ± 0.005 would yield a better fit to observational data (the WMAP value of r\ and 
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measured 2 H, 4 He, and 7 Li abundances) for big bang nucleosynthesis. This would correspond 
to an increase in the quark mass of SX q /X q = 0.014 ± 0.004 (here X q = m q / Aqcd)- 



Dent, Stern, and Wetterich 



111 ] calculated the sensitivity of BBN abundances for 2 H, 4 He 



and 7 Li to the variation of binding energies of 2,3 H, 3,4 He, 6,7 Li and 7 Be in a linear approx- 
imation. We use the response matrix values in their Table 1 for to B^ e and multiply 
by the (5mjv/mjv and K values of our Table IVl to estimate the sensitivity of BBN yields to 
variation of the quark mass. If we compare to the ratio of observation and theoretically 
predicted abundances given in their Appendix B, we obtain the following equations for 2 H, 
He and Li: 

l + 7.7x = 2,8± ° ,4 = 1.07 ±0.15 , (42) 
2.61 ±0.04 ' v ; 

0.249 ± 0.009 , , 

1 _ o.95x = = 1.005 ± 0.036 , 43 

0.2478 ±0.0002 ' v ' 

1 - 50s = L5 ± 0,5 = 0.33 ± 0.11 , (44) 
4.5 ±0.4 v ' 

where x = 5X q /X q . These equations yield 3 consistent values of x: 0.009 ± 0.019, —0.005 ± 
0.038 and 0.013 ± 0.002. The statistically weighted average of 5X q /X q = 0.013 ± 0.002 
is dominated by the 7 Li data. A more accurate calculation should take into account the 
effect of the 8 Be binding energy variation (which is not calculated in Ref. [ll(), the variation 
of the virtual 1 So(^p) level, and non-linear corrections in x which are important for 7 Li. 
Allowing for the theoretical uncertainties we should understand this BBN result as 5X q /X q = 
K ■ (0.013 ±0.002) where K ~ 1 and the expected accuracy in K is about a factor of 2. Note 
that here we neglect effects of the strange quark mass variation. A rough estimate of these 

Lin 

effects on BBN due to the deuteron binding energy variation was made in Refs. [9J, I10| . 



V. CONCLUSIONS 



We have argued that there are several reasons to question the spatial and temporal 
invariance of various fundamental "constants" of nature, such as the fine structure constant 
a or a comparable strong interaction parameter X q = m q j 'Aqcd- The search for evidence of 
such variations is ongoing in areas as diverse as quasar absorption spectra, the Oklo natural 
nuclear reactor, and big bang nucleosynthesis. 

In this work we have examined how nuclear binding energies depend on hadronic masses, 
including m^, ttia, rn-K, and a generic heavy meson my. We have done this by identifying the 
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mass-dependence in several realistic Hamiltonians — interactions that fit NN elastic scat- 
tering data and reproduce light nuclei binding energies reasonably well in quantum Monte 
Carlo calculations. By making small changes in the masses and re-evaluating the energy we 
have obtained the dimensionless derivatives of the energy with respect to variations in the 
hadronic masses. 

We have combined these results with a specific prediction from a Dyson-Schwinger equa- 
tion study of sigma terms in the light-quark hadrons for the hadronic mass-dependence on 
the quark mass m q . The pion mass changes most rapidly with changes in m q , so we find 
that both the one- and two-pion exchange parts of the NN interaction are very important 
for the consequent variations in nuclear binding. With our most complete model, the 1 So 
virtual bound state and deuteron vary in concert; if X q increases, they both become less 
bound, while if X q decreases, they both become more bound. (We note that this result is 
in disagreement with chiral perturbation results at the NLO level which have 1 S and 3 Si 
scattering varying antithetically.) The binding energies of A=3-8 nuclei behave in the same 
manner, all moving up or down together, with a sensitivity K in the range —1 to —1.5. 

Finally, we have folded these results with a study of the sensitivity of big bang nucleosyn- 
thesis to variations in nuclear binding. We find that a small increase in the quark mass of 
order 1% at the time of BBN is sufficient by itself to resolve existing discrepancies between 
theoretical and measured abundances of 2 H, 4 He, and 7 Li. 
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TABLE I: Ground state energies of light nuclei in MeV for the different Hamiltonians used in this 
work compared to experiment. 





1 S (np) 


2 H 


3 H 


3 He 


4 He 


5 He 


6 Li 


7 Li 


7 Be 


AV28 


0.0661 


-2.2250 
















AV14 


0.0663 


-2.2250 


-7.50 


-6.88 


-23.60 


-21.26 


-24.31 


-28.31 


-26.85 


AV18+UIX 


0.0665 


-2.2246 


-8.25 


-7.49 


-27.50 


-25.26 


-28.22 


-33.33 


-31.74 


Expt. 




-2.2246 


-8.48 


-7.72 


-28.30 


-27.41 


-31.99 


-39.24 


-37.60 



8 Be 

-40.26 
-48.50 
-56.50 



TABLE II: Dimensionless derivatives AS (ran ) = Sm H /m H °^ the ener gy f° r the 1 So( n P) virtual 
bound state e v and the deuteron Q for all three Hamiltonians. 



m H 




Ae„ 






AQ 




AV28 


AV14 


AV18 


AV28 


AV14 


AV18 


m N 


-88.1 


-32.6 


-33.4 


13.06 


8.63 


8.90 


mN + &N 




-91.2 


-121.2 




13.03 


17.82 


mA 


63.9 






-5.15 






5a 




68.1 


102.2 




-5.12 


-10.36 


m n (OPE) 


9.5 


-4.1 


-3.8 


-2.23 


-1.55 


-1.40 


m w (+TPE-s) 


24.4 


35.5 


53.0 


-3.63 


-4.02 


-6.70 


m w (+TPE-L) 








-4.02 


-4.31 


-6.74 


my 


-153.7 


-245.0 


-381.9 


20.88 


22.92 


41.74 
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TABLE III: A£(m H ) for the AV14 Hamiltonian. 





3 H 


3 He 


4 He 


5 He 


6 Li 


7 Li 


7 Be 


8 Be 


m N 


6.00 


6.44 


3.97 


4.58 


5.25 


5.60 


5.88 


5.10 


mN + 5n 


12.32 


13.17 


9.03 


10.38 


11.35 


12.74 


13.41 


11.71 




-7.35 


-7.82 


-5.89 


-6.74 


-7.10 


-8.31 


-8.76 


-7.69 


m w (OPE) 


-0.45 


-0.50 


-0.20 


-0.24 


-0.36 


-0.30 


-0.32 


-0.23 


m w (+TPE-s) 


-4.35 


-4.66 


-3.33 


-3.87 


-4.19 


-4.83 


-5.09 


-4.38 


m n (+TPE-L) 


-4.53 


-4.85 


-3.47 


-4.04 


-4.40 


-5.06 


-5.34 


-4.59 


my 


29.36 


31.30 


23.60 


27.09 


28.98 


33.72 


35.55 


30.98 


TABLE IV: A£(m H ) for the AV18+UIX Hamiltonian. 




3 H 


3 He 


4 Hp 


5 Hp 


6 Li 


7 Li 


7 Bp 


8 Bp 


m N 


6.07 


6.54 


3.99 


4.51 


5.12 


5.24 


5.49 


4.81 


mN + $N 


16.56 


17.73 


11.86 


13.31 


14.41 


15.53 


16.29 


14.36 


6a 


-12.20 


-13.02 - 


-9.16 - 


-10.24 


-10.80 


-11.96 - 


12.56 


-11.11 


m^ (OPE) 


-0.37 


-0.42 - 


-0.19 


-0.24 


-0.36 


-0.29 


-0.30 


-0.23 


m w (+TPE-s) 


-6.90 


-7.38 - 


-5.11 


-5.82 


-6.33 


-6.95 


-7.30 


-6.34 


(+TPE-L) 


-6.87 


-7.36 - 


-5.06 


-5.75 


-6.24 


-6.84 


-7.18 


-6.24 


(+TNI) 


-6.91 


-7.40 - 


-5.12 


-5.82 


-6.31 


-6.91 


-7.26 


-6.31 


my 


47.98 


51.23 


36.34 


40.87 


43.48 


48.11 


50.53 


44.40 


TABLE V: Dimensionless derivatives K 


- _ 5E/E 
~ &m q /m t 


- of the 


energy over light quark 


mass m q 


for the 


different Hamiltonians. 




















2 H 3 H 


3 He 


4 He 


5 He 


6 Li 7 Li 


7 Be 


8 Be 


AV28 


4.5 


-0.75 














AV14 


7.3 


-0.84 -0.89 


-0.96 


-0.69 


-0.81 


-0.89 -1.03 


-1.09 


-0.92 


AV18+UIX 


11.4 


-1.39 -1.44 


-1.55 


-1.08 


-1.24 


-1.36 -1.50 


-1.57 


-1.35 
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Fig. 1 (Flambaum & Wiringa) 




0.5 1 1.5 



m q /m q0 

FIG. 1: (Color online) Variation of two-nucleon energies with current-quark mass variation from 
DSE calculation: full calculation (solid lines) and with OPE modification only (dashed lines) 
for three different Argonne Hamiltonians. Virtual bound state energies are plotted as positive 
quantities; m q o is the physical current-quark mass. 
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Fig.2 (Flambaum & Wiringa) 
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Fig. 3 (Flambaum & Wiringa) 
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FIG. 3: (Color online) Variation of nuclear energies with current-quark mass variation from DSE 
calculation for AV18+UIX Hamiltonian. 
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